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The current understanding of mechanical feedback is reviewed by evaluating the standard, adi- 
_h ■ abatic model for shell formation and evolution. This model is relevant to phenomena ranging 

from individual stellar-wind bubbles to galactic superwinds, forming the basis for our under- 
standing of the multiphase ISM, IGM, and galactic evolutionary processes. Although significant 
discrepancies between the model and observation have been identified, to date there are none 
that require a fundamental revision. A variety of evidence, ranging over three orders of magni- 
tude in spatial scale, is broadly consistent with the standard model. This includes kinematics of 
t-H , individual objects, observations of hot gas, the size distribution of H I shells, and outflow rates 

from starburst galaxies. However, some of the most pressing issues relating to shell evolution 
are still outstanding and obstruct efforts to resolve key questions like the fate of the hot gas. 

8 ■ 1 Introduction 

£3 ■ 
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Mechanical feedback from massive stars is one of the principal drivers of evolutionary processes 
. in galaxies. Our understanding of the feedback process is grounded in the conventional model 
for the evolution of stellar wind- and supernova-driven superbubbles (e.g., Pikel'ner 1968; Weaver 
etal. 1977): hot (log T/K ~ 6 - 7), low-density (n ~ 0.01 cm" 3 ) gas is generated within a double- 
shock structure, and the pressure of this hot gas, chemically enriched by the stellar products, drives 
the growth of the thin, radiatively-cooled shell. In the standard model, energy loss from the hot 
gas is considered to be negligible, yielding simple analytic expressions for the shell radius R and 
expansion velocity v as a function of time t: 

R oc (L/n) 1 / 5 f 3 / 5 , 
v = R oc (L/n) 1/5 t- 2/5 . (1) 

For a stellar wind-driven bubble, the input mechanical luminosity L = 1/2 Mn^, where M and 
Voo are the stellar wind mass-loss rate and terminal velocity, respectively. For OB associations, 
supernovae quickly dominate over winds, in which case L = N*Es\/t e (e.g., Mac Low & McCray 
1988), where A* and E$\ are the total number of supernovae and supernova energy, respectively, 
and t e is the time over which the supernovae occur. 

The fate of this hot gas is crucial in understanding the phase balance and enrichment of the 
interstellar and intergalactic media. Is stellar feedback indeed the source of the diffuse, hot gas in 
the ISM? Do galactic superwinds from starbursts eject gas from galaxies? How does mechanical 
feedback affect the structure of the ISM and porosity, for example, to ionizing radiation? Do 
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superbubbles trigger renewed star formation? To answer these basic questions, we must first 
evaluate the relevance of the standard, adiabatic model for shell evolution, which is applied over 
scales ranging from single stellar wind bubbles to galactic superwinds. 

2 Single star bubbles 

Single star bubbles of isolated OB stars are woefully understudied. Oey & Massey (1994) examined 
two nebular examples in M33, and spectroscopically classified the late-type O stars. The inferred 
stellar masses and ages implied wind parameters that were consistent with the observed sizes and 
shell ages predicted by the adiabatic model. However, the parameters were loosely constrained 
and the shell kinematics were not observed. H I shells with radii of several tens of pc have been 
identified as wind-blown bubbles around a number of Galactic O and Of stars (Cappa & Benaglia 
1998; Benaglia & Cappa 1999). These are largely consistent with the standard model, and probe 
a specific subset of fairly evolved stars with old shells that have essentially stopped expanding. 

Studies of Wolf-Rayet ring nebulae suggest shells that are too small, equivalent to an overesti- 
mate in L/n by an order of magnitude (e.g., Treffers & Chu 1982; Garcfa-Segura Sz Mac Low 1995; 
Drissen et al. 1995). However, the progenitor star produces several wind phases, including both fast 
and slow winds, with strongly changing M. Their cumulative effect on shell morphology is even 
more complex and poorly-understood than the wind evolution itself, and it is therefore unsurpris- 
ing to find significant discrepancies between the predictions and observations of shell parameters. 
Hence it is desireable to study more of the simpler, single OB star bubbles. 

3 Superbubbles 

Superbubbles created by OB associations are more prominent than single-star bubbles, and thus 
have been studied more extensively. Soft X-ray emission has been detected from the interiors of 
many objects in projection, which is qualitatively consistent with the adiabatic superbubble model. 
Quantitatively, two classes of X-ray emission have been identified: objects with X-ray luminosity 
L x in excess of the standard model's prediction (Chu & Mac Low 1990; Wang & Helfand 1991), 
and objects that remain undetected in X-rays (Chu etal. 1995). The X-ray-bright objects are 
suggested to be overluminous because of SNR impacts on the shell walls. Upper limits on the X- 
ray-dim objects remain consistent with L x predicted by the adiabatic model. It will thus be of great 
interest to determine L x for these objects with XMM-Newton or Chandra. The existence of hot 
gas within superbubbles also implies an interface region with the cooler shells where intermediate 
temperatures and ions should be present. Chu et al. (1994) searched a number sightlines through 
LMC superbubbles and confirmed the existence of Civ and Siiv absorption in all cases. 

A stringent test of the adiabatic model is to compare the predicted and observed shell kine- 
matics where the input mechanical power and other parameters are well-constrained. We obtained 
spectroscopic classifications of the OB associations in eight LMC superbubbles (Oey & Massey 
1995; Oey 1996; Oey & Smedley 1998), thereby providing estimates of L in these young, wind- 
dominated objects. The predicted growth rate for the shells was higher than implied by their 
observed R and v, equivalent to an overestimate in L/n by an order of magnitude. However, even 
after adjusting L/n in the models, over half the objects still showed observed expansion velocities 
that were typically a factor of two higher than predicted for the given R. Similar discrepancies 
were reported for Galactic objects by Saken etal. (1992) and Brown etal. (1995). Figure [l] shows 
representative examples: the two objects in the left column show predicted R (solid curves) and 
v (dotted curves) in agreement with the observed values (horizontal segments); while the objects 
on the right show anomalously high observed v. All the models have L/n reduced by a factor of 
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Figure 1: Models for the evolution of LMC superbubbles, with L/n reduced by a factor of 10 from the 
empirical estimate. Solid, dotted, and dashed lines show the predicted R (pc), v (kms -1 ), and input 
10 log i/(10 35 ergs -1 ), respectively, as a function of t (Myr). The horizontal solid and dotted lines show 
the observed R and v, respectively, with the time constraints set by the observed stellar population. 



10. The objects with anomalously high v also each exhibit anomalously high L x and [Sn]/Ha, 
consistent with the suggestion by Chu & Mac Low (1990) and Wang & Helfand (1991) that SNR 
impacts have heated the shell. These would then also accelerate the shell walls, explaining the high 
observed v. The observed stellar mass function in these objects also support this scenario, implying 
the former existence of 1-2 higher-mass stars corresponding to SN progenitors. 

While the SNR impact hypothesis appears the most likely explanation for the high-velocity 
shells, a sudden drop in the ambient density can induce a "mini-blowout" with shell kinematics 
that can easily reproduce the anomalous observations (Oey & Smedley 1998; Mac Low et al. 1998; 
Franco & Silich 1999). Indeed, were it not for the high X-ray emission and additional evidence 
above, it would be impossible to distinguish between the two shell acceleration mechanisms from 
the kinematics alone. 

We therefore see that the ambient density distribution is critical in determining the shell evolu- 
tion. An underestimate in n could, for example, contribute to the growth rate discrepancy described 
above, that is seen in all the objects. To clarify the effect of the ambient gas distribution, Oey 
et al. (2001) mapped the H I distribution within a ~ 40' radius of three nebular LMC superbubbles 
at 30" resolution, using the Australia Telescope Compact Array. The results show neutral environ- 
ments that vary to an extreme: DEM L25 shows no H I shell component, but appears to be nestled 
in an H I hole; DEM L50 shows an H I shell with ~ 10 times more mass than the Ha component, but 
is otherwise in a large (1.4 kpc) Hi void; and DEM L301 is in a non-descript neutral environment 
with no apparent correspondence to the optical shell. Thus, despite the similar appearance of the 
Ha nebulae, the H I distributions could not be more different. It is therefore extremely difficult to 
infer properties of ambient density distribution without direct, empirical determinations. 
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Another vital parameter for shell evolution is the ambient pressure, Pq, which determines 
whether and when the superbubble growth becomes pressure-confined. While Po is usually rel- 
atively unimportant in young, high-pressure superbubbles such as the nebular objects modeled 
above, it is of vital importance in the mid- to late-stage evolution. It may also be relevant in 
high-pressure, ionized environments like dense star- forming regions (e.g., Garcfa-Segura & Franco 
1996). The value of Po ultimately determines the final size of the objects, and conditions relative 
to blowout. The uniformity and distribution of Po i n the multiphase ISM is therefore especially 
relevant to the global effect of superbubbles on the ISM (see below) . 

Finally, if the hot gas within superbubbles does not somehow blow out and merge into the hot, 
ionized medium (HIM), it is likely that the objects will cool and depart from energy conservation. 
Indeed, cooling of the hot interior has long been one of the principal questions for superbubble 
evolution and the fate of the hot gas. Thermal conduction at the interface between the cool 
shell wall and hot gas should cause a high rate of mass-loading into the interior, by evaporation. 
The evaporated shell material strongly dominates the mass in the hot region, which could be 
further supplemented by evaporation and ablation from small clouds overrun by the expanding 
shocks (e.g., Cowie & McKee 1977; McKee etal. 1984; Arthur & Henney 1996). If the interior 
density is sufficiently increased, radiative cooling will dominate, and the shells will no longer grow 
adiabatically. In addition, Silich etal. (2001) point out the importance of enhanced metallicity 
in the superbubble interiors, caused by the stellar and SN yields. Preliminary investigation for 
individual objects by Oey & Silich (2001, in preparation) shows enhancement in L x by almost an 
order of magnitude for low-metallicity [Z = 0.05 Zq) objects. The elemental yields therefore may 
also significantly enhance the cooling rate, facilitating the transition from adiabatic to momentum- 
conserving evolution. 

4 Global mechanical feedback and superwinds 

Quantifying the statistical properties of superbubbles in galaxies allows us to test the evolutionary 
model and also gain insight on the properties of the ISM. In particular, the so-called interstellar 
porosity (see below) is a conventional parameter for evaluating the relative importance of the HIM 
in the multiphase ISM. 

4.1 Superbubble populations 

Oey & Clarke (1997) derived expressions for the differential size distribution N(R) dR of super- 
bubbles in a uniform ISM, using the analytic expressions for adiabatic evolution (equation ||). We 
considered a power-law mechanical luminosity function for the parent OB associations, 

(p(L) dL = AL-P dL , (2) 

with [3 = 2, based on the typical H n region luminosity function (H n LF; e.g., Kennicutt et al. 1989; 
Oey & Clarke 1998a). The superbubble growth was also taken to be pressure-confined when the 
interior pressure P, = Po. Star formation is assumed to be coeval within each OB association, with 
SNe therefore exploding over a period t e = 40 Myr, the lifetime of the lowest-mass SN progenitors. 
These simple assumptions yield two evolutionary paths for the superbubbles: objects whose growth 
stalls at radius Pf due to pressure confinement at time tf oc (n/Po) 1 / 2 Pf; and objects whose growth 
continues until t e . The characteristic timescale t e therefore determines the characteristic parameters 
R e and L e , which correspond to conditions for which an object stalls at t e . 

The resultant size distribution for a constant superbubble creation rate tp, and single- valued 
4>{L) is: 

N gTOW (R) oc P 2 / 3 . (3) 
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This is determined exclusively by growing objects, since the stall radius R{ in this case is the same 
for all the shells, and represents the maximum R in the distribution. 

If all objects are created in a single burst, but now with a power- law <j)(L) (equation |2|), then 
both growing and stalled objects contribute to the size distribution: 



N. 



grow(-R) OC R 



N staU (R) oc fr-WxR 



4-5/3 
1-2/3 



oc R~ 



(4) 



with (3 = 2. The stalled objects dominate the total N(R) for R < Rf(t e ), beyond which the growing 
objects dominate, falling off with a much steeper independence. 

For constant i/j and power-law <p(L), the stalled objects again dominate, with a total N(R): 



rl-/3 / r> 

N(R) = A^ 

K P \-tte. 



1-2/3 



2{t e +t s ) + 



9-6/? 

-2 + 3/r 



te. — 



R 

Re 



(5) 



where t s is an arbitrary period after t e during the which the superbubbles survive before disap- 
pearing. Here again, the stalled objects dominate, effectively yielding N(R) oc R l ~~ 2/3 oc R~ 3 . The 
maximum in N(R) is at Rf (L m in), corresponding to single-progenitor SNRs. At R < R{(L m i n ), 
growing, high-L objects dominate with N(R) oc R 2 ^ 3 (equation ||). Oey & Clarke (1998b) show 
schematic representations of N(R) for the above three combinations of ip and <f>(L). 

We compared this result to the H I shell catalog for the Small Magellanic Cloud (SMC) compiled 
by Staveley-Smith etal. (1997). This is by far the most complete sample of Hi shells obtained for 
any galaxy, as evidenced by the fact that the relative number counts of H n regions and H I shells 
are in excellent agreement with their relative life expectancies. For shells having R > 100 pc, the 
fitted power-law slope a = 1 — 2(3 is 2.7 ± 0.6, in excellent agreement with the general prediction 
of a = 3. Using (3 = 1.9 from the observed Hn LF for the SMC (Kennicutt etal. 1989) yields a 
specific SMC prediction of a = 2.8 ± 0.4, still in remarkable agreement with the observed value 
(Oey & Clarke 1997). 

We note that different models for ISM structure yield different predictions for N{R). For 
example, Stanimirovic et al. (1999) suggest a possible fractal structure for the neutral ISM. From 
the same H I dataset of the SMC, they find a fractal dimension implying a size distribution for H I 
holes of a = 3.5. Within the uncertainties, this prediction cannot be empirically differentiated from 
our model with a = 3; but it is worth noting that the predictions are in fact intrinsically different. 

However, the superbubble size distribution presently is not a sensitive test in determining 
whether or not the objects evolve adiabatically. If all the internal energy is radiated away, the 
objects are predicted to follow the momentum-conserving law given by Steigman etal. (1975): 



R oc {L/n Voo ) l/A t 1 ' 2 



(6) 



The stall radius R{ in this case is only 1.3 times larger than for the adiabatic model. Furthermore, 
the size distribution follows the same law, N(R) oc R 1 - 213 (Oey & Clarke 1997). The observations 
of hot gas are therefore vital confirmation that the adiabatic model applies to a significant fraction 
of superbubbles. 

It is also possible to derive the distribution of expansion velocities N(v) dv. This is naturally 
determined only by the growing objects (Oey & Clarke 1998b): 



iVgrow(^) OC V- 7 ' 2 , (3 > U 



(7) 



We again compared this prediction to the SMC H I shell catalog, and found a fitted power-law slope 
of 2.9 ± 1.4, which is also in agreement with the prediction. Thus, in spite of the crude assumptions 
in deriving the shell size and velocity distributions, the data suggest that the neutral ISM in the 
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Table 1: Porosities in the Local Group 



Galaxy 


Q 


M31 


0.03 


M33 


0.3 


LMC 


1.0 


SMC 


0.3 


IC 10 


23 


LG dwarfs 


0.01 - 0.2 


Milky Way (H n LF) 


0.2 


Milky Way (SN rate) 


-1 



SMC is fully consistent with superbubble activity dominating the structure. Although most other 
available H I shell catalogs are highly incomplete, preliminary results for a few other galaxies also 
show agreement with our model for the size distribution (Kim etal. 1999; Mashchenko etal. 1999; 
Oey & Clarke 1997). 

4.2 ISM porosity and galactic superwinds 

It is straightforward to use the analytic expression for N(R) (equation ||) to derive the interstellar 
porosity parameter Q, which is given by Oey & Clarke (1997) for 2- and 3-dimensional cases. Q 
is the ratio of: total area or volume occupied by superbubbles, to the total area or volume of the 
galaxy. Thus it is essentially the filling factor of hot gas, assuming hot gas is contained within 
all of the superbubbles. Determining Q has been the conventional way to evaluate the relative 
importance of the HIM (e.g., Cox & Smith 1974; McKee k Ostriker 1977; Heiles 1990). Values 
near unity indicate the HIM dominates the multiphase ISM, and values significantly > 1 imply an 
outflow, with the galaxy generating more hot gas than it can contain. 
We can write Q in terms of a galaxy's star- formation rate (SFR), ^f: 

vfr(M yr- 1 ) 1 

Q " 16 K(kpc3) °<J%' (8) 

for (3 = 2, a Salpeter (1955) IMF for stellar masses 0.1 < m < 100 M , and P /k = 9500. R g and 
h are the galaxy's gaseous radius and scale height, respectively. We caution that Q depends on 
ambient interstellar parameters, for example, Pq 1 as indicated. 

Table 1 shows the interstellar porosities for all the galaxies in the Local Group. Dwarf galaxy 
parameters are taken from Mateo (1998), and SFR's for M31, M33, and the Magellanic Clouds 
from their Hn LF's (Kennicutt etal. 1989). In general, we see that almost all the galaxies show 
low values with Q < 0.5, implying that the HIM is not the dominant ISM component in these 
cases. The LMC is an exception, showing Q ~ 1, so the HIM appears likely to dominate in this 
galaxy. The Milky Way results historically have been ambiguous and controversial (e.g., McKee & 
Ostriker 1977; Slavin & Cox 1993), and we obtain different answers from different methods. Taking 
the Galaxy's SFR from the Hn LF of McKee & Williams (1997) and Smith & Kennicutt (1989) 
yields a low Q ~ 0.2; but using the SN rate of van den Bergh & Tammann (1991) yields Q 1. 

However, by far the highest porosity is found for IC 10, the one starburst galaxy in the Local 
Group. IC 10 shows a Q value that is over an order of magnitude larger than for any other galaxy 
examined, therefore fulfilling the breakout criterion for a galactic superwind outflow. It is possible 
to estimate the mass-loss rate in this outflow M out , assuming that the material is largely evaporated 
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from shell walls by thermal conduction. Mac Low & McCray (1988) have estimated the evaporation 
rate for an individual superbubble: 

M(R) ~ lQ-K f im ll /25k CT 5 ^ 2 R , (9) 

where /imn is the mean mass per nucleon and C is a conductivity coefficient. Integrating this over 
the superbubble size distribution yields M out : 

M out =* / M(R) N(R) dR . (10) 

The estimate of M ou t is necessarily crude, for example, we rather arbitrarily take R max = R e , with 
the caveat that it may be limited by R g . 

The maximum M out is obtained for the largest evaporating surface area of superbubbles, which 
occurs for Q ~ 1. Porosities of Q > 1 imply strong interaction and destruction of the shells, or 
blowouts, probably reducing the available surface area. For IC 10, the derived normalization for 
N(R) yielding Q = 1 corresponds to M out ~9x 10~ 3 M Q yr _1 . For comparison, the SFR required 
for Q = 1 (equation |8|) is \P ~ 3 x 10~ 2 M Q yr _1 , which is only about 3 times more than M out ! 
Although the actual SFR for IC 10 is much larger (7 x 10 _1 M yr _1 ), yielding the higher original 
Q value, this rough calculation suggests that M out from galactic winds can be of similar magnitude 
as the parent starburst SFR. Indeed, an extensive absorption-line study of local starburst galaxies 
by Heckman et al. (2000) also concluded that empirically, the outflow and star-formation rates have 
the same order of magnitude for their sample. 

Finally, we note the extensive body of numerical work on superbubbles and blowout conditions. 
These are too numerous to discuss here, but we refer the reader to Mac Low (1999) and Strickland 
&i Stevens (2000), who provide overviews of this field. The details of the numerical predictions are 
presently difficult to confirm empirically, but observations with XMM-Newton and Chandra will be 
especially helpful to constrain the dominant processes. 



5 Summary 

We have seen that observations of mechanical feedback ranging from individual stellar wind bub- 
bles to galactic superwinds are all roughly consistent with the conventional adiabatic model for 
bubble evolution. The empirical confirmation of hot, coronal gas supports the relevance of the 
adiabatic evolution for at least an important proportion of shell phenomena. Although kinematic 
discrepancies are often found for individual objects, these anomalies have plausible explanations, 
and presently there are no cases that suggest a need for any major revision of the conventional un- 
derstanding. It is a remarkable strength that the adiabatic model apparently applies to phenomena 
over a range of size scales covering at least three orders of magnitude, from pc to kpc. 

However, characterizing the dominant parameters and their effects on the shell evolution is 
still highly problematic. For example, critical ambient ISM conditions like the density, pressure, 
and ionization distributions remain elusive. The mechanisms and conditions for cooling of the 
interior energy need to be identified, and the energy budgets reliably determined. Perhaps the most 
fundamental question is the fate of the hot gas generated within the superbubbles: Does it escape 
and thereby constitute the HIM? Does it indeed escape from starburst galaxies, and from their 
gravitational potentials? These issues have vital consequences for almost all galactic evolutionary 
processes, and our understanding depends on further clarifying the bubble and superbubble process. 

MSO gratefully acknowledges travel support from the Graduiertenkolleg and STScI Institute 
Fellowship. Some material in §4.2 was derived in 1998 at the Aspen Center for Physics. 
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